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The derivation of the effective theory for the phase degrees of freedom in a superconductor is 
still, to some extent, an open issue. It is commonly assumed that the classical XY model and its 
quantum generalizations can be exploited as effective phase-only models. In the quantum regime, 
however, this assumption leads to spurious results, such as the violation of the Galilean invariance in 
the continuum model. Starting from a general microscopic model, in this paper we explicitly derive 
the effective low-energy theory for the phase, up to fourth-order terms. This expansion allows us 
to properly take into account dynamic effects beyond the Gaussian level, both in the continuum 
and in the lattice model. After evaluating the one-loop correction to the superfluid density we 
critically discuss the qualitative and quantitative differences between the results obtained within 
the quantum XY model and within the correct low-energy theory, both in the case of s-wave and d- 
wave symmetry of the superconducting order parameter. Specifically, we find dynamic anharmonic 
vertices, which are absent in the quantum XY model, and are crucial to restore Galilean invariance 
in the continuum model. As far as the more realistic lattice model is concerned, in the weak-to- 
intermediate-coupling regime we find that the phase-fiuctuation effects are quantitatively reduced 
with respect to the XY model. On the other hand, in the strong-coupling regime we show that the 
correspondence between the microscopically derived action and the quantum XY model is recovered, 
except for the low-density regime. 

I. INTRODUCTION 

The description of the low-energy dynamics of the phase 9 of the superconducting order parameter A = |A|e*^ has 
recently attracted a renewed interest in connection to the phenomenology of high- Tc superconducting cuprates. These 
materials are characterized by a strongly anisotropic (quasi-two-dimensional) crystal structure and by anomalous 
normal-state properties, which can be interpreted in terms of preformed Cooper pairs (non-zero amplitude |A|) 
without phase coherence (vanishing superfluid density Ps)^ Although in ordinary (weak-coupling) superconductors 
both |A| and ps vanish at the critical temperature Tc, due to thermally excited quasiparticles, it is possible that ps 
vanishes when |A| is finite, due to phase fluctuations, ff such is the case, phase fluctuations should play a role more 
relevant than expected in ordinary superconductors, both near Tc and deeply in the superconducting state. 

A crucial point in the analyses of phase-fluctuation effects is clearly the choice of the effective model, the most com- 
mon being the classical XY model, as the paradigm for the universality class which is relevant to superconductivity!^ 
The XY model, which has also been extensively used in contexts such as systems of resistively shunted Josephson 
junctions or granular superconductors^^ assumes a phase field Oi defined on a coarse-grained lattice with a lattice con- 
stant determined by the coherence length which sets the length scale above which the fluctuations of the amplitude 
|A| become uncorrelated, and |A| may be assumed to be fixed. Thus, at distances larger than phase fluctuations 
only are relevant, and can be described in terms of a Josephson-like interaction between neighboring sites < i,j >, 
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Here 6*^ = 6i ~ 9j ~ |V6'|^o is the phase difference between nearest-neighboring sites, and is the bare coupling 
constant. Eq. is written for a two-dimensional system (which is the case relevant for cuprate superconductors), 
fn d = 2, is an energy which measures the superfluid stiffness, i.e., the energetic cost to produce phase variations 
in the system. In the continuum (Galilean invariant) model Dq = ps{T = 0)/m — p/m, where m is the electron mass, 
p is the electron density, and we set h — 1. In the case of a generic d-dimcnsional system, we still refer to Dq as to 
the stiffness, even though the corresponding energy scale is Z?oCo~^- 

Within the XY model, one can determine the depletion of the superfluid stiffness due to the collective-mode 
excitations, i.e., to the anharmonic terms in the gradient expansion of Eq. (Q. At low temperature, the leading 
correction to the bare stiffness Dq comes from the fourth-order term, and within perturbation theory with respect to 
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the harmonic term one gets, e.g. in d = 2, 
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where Dq = D{T = 0), since ((V6')^) = 4r/^gZ?o when evaluated at Gaussian level. This result led to the proposal 
that phase fluctuations, rather than the d wave quasiparticles excitations characteristic of the superconductivity in the 
cuprates, are responsible for the linear thermal depletion of the superfluid density^ which is experimentally observed 
in these materials, down to very low temperatures (1 — 5 K)4i However, one would expect that quantum effects are 
relevant at such low temperatures, and, even though the XY model may be reasonably adopted in the classical regime, 
the investigation of the quantum regime is much more involved. 

It has been proposed^ that the quantum effects can be partially included in the so-called quantum XY model, 
by deriving the Gaussian phase-only effective action from a microscopic BCS model, while obtaining the anharmonic 
terms through the expansion the cos^y term of Eq. J^l- As a consequence, the interaction terms in the phase, given 
by powers of 9ij ~ Co|V0|, are purely classical, whereas the Gaussian propagator is evaluated in the quantum regime. 
However, this approach turns out to be unsatisfactory within many respects. 

Indeed, since in the quantum case ((V6')^) ^ aX T — Q (see below), Eq. (O leads to a finite correction to Dq 
even at T = 0, which, in turn implies Ps(0) 7^ p in the continuum model, thus explicitly violating Galilean invariance. 
Therefore, the description of the phase fluctuations by means of the quantum XY model misses some important effect 
within the continuum model. A second issue, which is particularly relevant in connection to the description of high- 
Tc superconductors, concerns the limit of strong pairing interaction. Indeed, while at weak coupling the amplitude 
fluctuations can be safely neglected in deriving the phase-only theory, as the interaction increases, the fluctuations of 
the modulus of the order parameter become intimately connected with the density fluctuationssiA Since the phase and 
the density are conjugate variables, the description of phase fluctuations requires a proper treatment of the amplitude 
fluctuations in the strong-coupling regime. 

In this paper we present a systematic comparison between the phase-only model derived microscopically and the 
quantum and classical XY models. To this extent, after describing the formal steps to derive the effective action in 
various cases (continuum vs lattice model, neutral vs charged system), we evaluate the one- loop corrections to _Do 
coming from the non-Gaussian (anharmonic) terms. 

We find that quantum effects generate dynamical interaction terms, absent in the quantum XY model, whose role is 
crucial, e.g., in preserving the Galilean invariance in the continuum model, and in reproducing the expected mapping 
of the lattice model onto the XY model in the strong-coupling regime. Some of the results on quantum corrections 
in the weak-coupling limit where already presented in a short paper and we provide here a detailed derivation, since 
the issue, which involves many subtleties, has been only partially addressed in the literaturei^4S4LS424i2iiiiiSiiiiiii& We 
mainly focus on the quantum (T = 0) corrections to the superfluid density, and we comment on the possible extension 
to the classical regime. Indeed, besides the restoration of the Galilean invariance in the continuum system, the 
microscopic derivation of a phase-only theory, developed in this paper, allows to demonstrate clearly that, in general, 
the XY model leads to an overestimation of the depletion of the superfluid 5D stiffness, even when the realistic effects 
such as the Coulomb screening are taken into account, with the only exception of the extreme strong-coupling regime 
where the XY and the microscopic estimate of 5D tend eventually to coincide. For sake of clarity we summarize 
below our final results, with reference to the corresponding Sections where they are derived, as a guideline for the 
reader. We report in tab. I and II (respectively for the continuum and the lattice models) the estimations for the 
ratio between the corrections to the superfluid density, derived within the microscopic models {5D) and within the 
XY model {5Dxy) in all the cases considered in this work. 



Continuum models 





Quantum regime (T = 0) 


Classical regime (high T) 




Neutral system 


Charged system 


Neutral system 


Charged system 


5D/5DxY 




(IV B) 




(IV B) 


(IV C) 


ep/ec ■ l/(fcF^o)''+' 
(IV C) 



TABLE I: Ratio between the superfluid-stifTness correction derived within continuum microscopic models and the XY model: 
the classical and the quantum case are explicitly considered. Here kp and ef are the Fermi momentum and the Fermi energy 
respectively, whereas ec is a typical Coulomb energy scale (cfr. Sec. IV C). Note that, quite generally, l/(fcFCo) << 1, i.e. the 
XY model tends to overestimate 5D). 
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Lattice models (neutral and charged system, quantum and classical regime) 





Weak Coupling 


Strong Coupling 




^ l/(fcFeo)^ 

(V) 


~ 1 
(VI) 



TABLE II: Ratio between the superfluid-stiffness correction derived within lattice microscopic models and the XY model: the 
classical and the quantum case are in this case qualitatively equivalent. 



The specific is sue of the temperature dependence of the superfluid density in high- Tc superconductors has been 
analyzed in Ref. |5l6| . In general, according to Eq. (Q) , one would expect that Coulomb effects, lifting the phase mode 
to the plasma frequency (see Sec. II B), completely suppress the contribution of phase fluctuations to D{T). However, 
in Ref. Q it has been shown that the d-wave symmetry of the order parameter plays a crucial role in determining the 
behavior of the phase mode in the presence of both long-range Coulomb interactions and dissipation. Even though 
here we ignore dissipative effects, for the sake of completeness we still compare s-wave and d-wave superconductors. 
Once again, the formalism of the effective action makes this comparison, to some extent, particularly simple, and of 
transparent physical interpretation. 

The plan of the paper is the following. In Sec. Ill Al we derive the Gaussian phase-only action starting from the 
continuum model, in the absence of long-range Coulomb forces. The case of the charged system is dealt with in Sec. 
IIIBI Although most of the results presented in Sec. ^Ja-re standard, our scope here is to introduce the formalism 
and a classification scheme for the interaction vertices, which turns out to be particularly useful in dealing with 
non-Gaussian terms. In Sec. IIIII we apply the same procedure to a lattice model, which does not appear as a trivial 
extension of the continuum model. In Sec. IIVI we calculate the anharmonic terms, providing a classification of the 
diagrams (Sec. II V A|) . and determine the one- loop correction to the stiffness within the continuum model, both at 
T — (Sec. IIVB|) and in the classical limit (Sec. IIV C|l . The results are then compared with the results of the XY 
model. In Sec. El the same analysis is carried out in the lattice case. The strong-coupling regime of the lattice model 
is discussed in Sec. IVII and the results are summarized in Sec. IVIII The details on the derivation of the phase-only 
action are reported in App. ^ whereas in App. 0we discuss the connection between the coefficients of the Gaussian 
action and the gauge-invariant electromagnetic response functions. 

II. THE EFFECTIVE ACTION FOR THE CONTINUUM MODEL 
A. The neutral system 

To introduce the general formalism adopted in this paper, and for the sake of completeness, we shortly discuss the 
standard formal steps to derive the phase-only action within the continuum BCS model H = Hq + iJ/, with 

^0 = -E/ rfr^c+(r)V2c.(r), (3) 
Hi = ^W"'(k')Ck+a|C+k+a4C_k' + axCk. + |t. (4) 

k,k',q 

Most of our results are obtained for generic spatial dimension d of the system, unless explicitly indicated. In Eqs. 
©-(01) , Cct^'' is the annihilation (creation) operator of an electron of spin cr, fl is the volume, and C/ > is the pairing 
interaction strength. Hereafter we assume H = = I- The factor w(k) controls the symmetry of the Cooper-pair 
wave function, e.g., w(k) = 1 for s-wave superconductors, and u'(k) = cos2(/), with (j) = aian^ky / k^) , for d-wave 
superconductors, in the continuum case and for d = 2. Note that, in the case of a d-wave symmetry, four nodes 
are present in the quasiparticle excitation gap, and the existence of gapless excitations is relevant in determining the 
low-temperature thermodynamics of the system. 

We first discuss the case of a neutral system, while the Coulomb interaction will be introduced in Sec. Ill Bl to 
deal with the charged case. As customary, the microscopic effective model for the collective modes is derived by 
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considering the action corresponding to the Hamiltonian (O-Q, within the finite-temperature Matsubara formahsm, 

S„^^cro ^So + Sl= /'^ dr i V (t) [Or + ^k]cka (T)dr + Hj (t) l , (5) 



where r is the imaginary time, [3 — \/T, and Ck^A* is the band dispersion with respect to the chemical potential 
/J,, with = k^/2m in the continuum case. To obtain the effective action in terms of the order-parameter collective 
degrees of freedom, the interaction is decoupled in the particle- particle channel by means of the Hubbard-Stratonovich 
transformation, introducing the auxiliary complex field A(x,t), as explicitly reported in App. ^ It is then possible 
to make the dependence on the phase of A(q, t) explicit in the action, by performing the gauge transformation Eq. 
(|A2|) on the fermionic fields. Indeed the dependence on the phase is eliminated from S*/, and made explicit in the free 
part of the action 5*0, which now reads 

^0 = 5-0+ ( dx (j)+ix)t(t){x), (6) 



with 
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where fi are the Pauli matrices and the operator V= (V — V) acts on the fermionic Nambu spinor (l){x) defined as 
the column vector (c|(a;), c^{x)). It is worth noting that the field 9 appears in the action only through its time and 
spatial derivatives, as it is expected for a Goldstone (massless) field. Moreover we point out that, while in the s-wave 
case the transformation (|A2|) completely eliminates the phase 9 from Sj, in the d-wave case a residual dependence 
survives. Nonetheless, in App. ^ where we report some more detail, we show that this residual dependence can be 
safely neglected in the hydrodynamic limit, which is the relevant regime when discussing the low-energy properties of 
the collective mode. Thus one can now integrate out the fermions, which appear quadratically in So + Sj, leading 
to the effective action (jA7|l . written in terms of the collective variables only. This standard procedure allows one 
to analyze the starting problem ©-Q from a different point of view. Instead of perturbatively expanding the 
Hamiltonian in powers of the coupling U, the effective action (|A7|) is expanded by assuming that the fluctuations 
of the fields around their saddle-point values are small, and that the variations of the phase in space and time are 
slow. Moreover, if one is interested in the dynamics of the phase at low temperature, the fluctuations of the modulus 
|A| can be neglected. As far as the specific problem of high-Tc superconductors is concerned, this approximation 
seems particularly appropriate for underdoped cuprates, which exhibit an anomalously large binding energy when 
compared to the critical temperature^ On the other hand, one should expect that modulus fluctuations can still be 
relevant at distances smaller than the typical length scale over which |A| fluctuates. It is then natural to introduce a 
spatial cut-off for phase fluctuations, of the order of the coherence length ^o- At distances greater than the phase 
fluctuations are the only relevant degrees of freedom, and a description in terms of the phase-only action becomes 
meaningful. Accordingly, the cut-off for the phase-only action is deflned as the characteristic length scale of the 
spatial decay of the correlation function for |A|. This point is delicate, since the naive identification of the coherence 
length ^0 with the Cooper-pair size S.pair is meaningless in c?-wave superconductor, where Cpair diverges at T = 0, due 
to the presence of the gapless quasiparticle excitations mentioned above, as it is discussed in Ref. 15]. 

Starting from Eq. I|A7|) of App. ^ with |A| = |A|(0) + (5|A|(q), and neglecting the fluctuations of the modulus 
(5|A|(q), the sum U^^ J2q I^K?)^ reduces to its g = value. As a consequence, the self-energy HA8|I does no longer 
depend on |A|. The action (|A7|) is then decomposed as Sneutrai — Smf + Seff{9), where Smf = I^P(0)/C^ — 
Tr [in 5,^^] is the mean-field BCS action, and Seff is the phase-only 

OO -. oo 

S.ff{9) = Tr ^ ^(E^o)^ = E ^S' (8) 

N=l N=l 

which takes into account phase fluctuations around the BCS saddle point. Here Qq is the BCS matrix defined in Eq. 
[9I) . and the trace Tr is explicitly given in fc-space by 

Sijj = ^ E tr[5^o(fci)S(fci - ka) X go(k2)t{k2 - fca) x Goik^MkN - ki)], 

All , . . . ,fciv 
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where tr is the trace in the Nambu space only and £(fc) is the Fourier transform of the expression Since the 
self-energy Q contains powers of the field 9 up to the second, the expansion ||SJ) is not directly an expansion in powers 
of 9 and the terms with the same power of 9 must be recollected, leading to the expansion 

oc oo 

Seff{0)=Y.^'^^J2 E Mqi,---,qk-l)0{qi)---9{qk-,)9{-q, qk-i), (9) 

k=2 k=2qi,---,qk-l 

which is completely equivalent to Eq. (jHJ, and includes both the Gaussian (fc = 2) and the anharmonic {k > 2) terms. 

Each term in the effective action (O is given by a closed fermionic loop Ak with k incoming lines corresponding to 
fluctuating 9 fields. In forming the fermionic loops, according to Eq. |(7J), we have three kinds of fcrmion-^? vertices, 
depicted in Fig. ^ 

(a) The vertex with a single incoming dotted line corresponds to the insertion of the time derivative of the phase, 
i.e., the term t]^ = {idr9/2)T3 of Eq. 10. 

(b) The vertex with two incoming dashed lines corresponds to the insertion of the gradient squared of the phase, 
i.e., the term t% = [{V 9f / 8m]T3 of Eq. (0. 

(c) The vertex with a single incoming dashed line corresponds to the insertion of the spatial derivative of the 
phase associated to the spatial derivative of the incoming and outgoing Green functions, i.e., the term T,]^ — 

i{[^9{x)- V]/4m}fo of Eq. ©. 

\ / I 

\ / I 
\ / I 

■ \ / I 

■ \ / I 

\ / I 



(a) (&) (c) 



FIG. 1: Feynman diagrams which correspond to the various insertions of the phase field according to Eq. (|7J. Dotted and 
dashed lines represent, respectively, the insertion of a time derivative or a gradient of the phase filed 6 in the fermionic lines 
(depicted with solid lines) 

The terms generated in the effective action by the self-energy = -I- E^, associated to the fa matrix, are 
essentially the same which are met when deriving the phase-only action in a bosonic superfluid system.^'' We therefore 
refer to these terms as "bosonic" contributions. The T,p term generates instead contributions which are strictly related 
to the existence of quasiparticle excitations. Therefore, we refer to them as "fermionic" contributions. Clearly, the 
expansion JSJ also generates "mixed" contributions, obtained by combinations of tq and terms of the self-energy 
E. Such a distinction will be relevant in the following, while discussing the derivation of the anharmonic terms of the 
phase-only action. 

At Gaussian level (fc — 2) the effective action describes the collective mode of the system. Beyond the Gaussian level 
(fc > 2), interaction terms in the field 9 appear, and we obtain the quantum generalization of the anharmonic terms of 
the XY model, Eq. (0. As in the classical case, the anharmonic terms are responsible for the reduction of the stiffness 
due to phase fluctuations. Before discussing this issue, we analyze the structure of the Gaussian phase propagator in 
the quantum regime. The Gaussian effective action is derived with long but straightforward calculations, and reads 

SZutrai = IJ2 l^'m^pi^) " <la^bAfj{q) + 2zr!„q, A" ^ (g)] \9{q)\\ (10) 
g 

The coefficients App, Apj and Ajj, which are depicted in Fig. Inland whose explicit expressions are reported in App. 
IXI represent respectively the density-density, density-current, and current-current bubbles, evaluated with the BCS 
Green functions. Since the BCS approximation explicitly breaks the gauge invariance^ the A functions do not 
coincide with the physical correlation functions K^^, which must obey the gauge invariance. However, as it is shown 
explicitly in App. ^ (where the general relationship between the A and the K bubbles is obtained), it is possible to 



6 




FIG. 2: Coefficients of the Gaussian action llUl . The first diagram comes from the = 1 term of Eq. Q, while the others 
come from the N — 2 term of Eq. Q. The correspondence between different linestyles and the various insertions is the same 
as in Fig. Q Notice that the full current-current correlation function Ajj — —pMp/mS"^ + AJ^ is given by the sum of the 
first and last diagram. 



relate the coefficients of Gaussian effective action with measurable quantities such as the bare compressibihty k(T) 
and the stiffness D{T) — ps{T)/m of the system, in the hydrodynamic limit^ Indeed, in this regime, we obtain the 
expression 

Shydr = \Y. [«^o(T)f^L + D{TW] d{q)e{-q). (11) 



The Bogoljubov sound mode of the neutral system, with velocity Cg = ^ D{T) / kq{T), is thus recovered. It is worth 
noting that, since the phase and the density are conjugate variables, the collective mode which appears in the phase- 
field propagator is the same as the density mode (see, e.g., Ref. (20) ), as it can be explicitly checked by deriving the 
physical density-density correlation function Kpp according to Eq. I|B1|) . 

By comparing Eq. (|ll|l with the classical XY model, one can recognize the same term ps(V9)^ which appears in 
Eq. (0). However, the quantum Gaussian propagator for phase fluctuations contains also the fi^ term of Eq. l(TT|l 
which takes into account the dynamic effects. Moreover, according to the results of App. iBland to Eqs. (jA13|l - ljA14|l . 
the temperature dependence of the bare stiffness DiT) which appears in Eq. 1111 is entirely due to the quasiparticle 
excitations. Indeed, the coefficient A^^ in Eqs. l)A13p - (|A14|l describes the normal-fluid component, whose contribution 
increases as the temperature increases 

k 

In the absence of dissipation the coefficients A are not analytic aX q — Q and finite temperature, since in this case the 
static (r^m = 0, q ^ 0) and the dynamic limit (q = 0, Vim 0) are different. This issue has been recently addressed 
in Ref. jl3| • Since we are interested in the T = correction to the superfluid density we shall not discuss further the 
consequences of such non analyticity, which would disappear in any case in the presence of dissipation. 



B. The charged system 

When dealing with a charged system, the long-range Coulomb interaction between the fermions, 

Hcoul^^ Xl^(^)^J+q,<T^^'-q.<T'Ck',a'Ck,<T, (13) 

must be included in the Hamiltonian ©-lEl)- V{q) is the Fourier transform of the three-dimensional Coulomb 
repulsion V{r) = e^/es?", projected onto a d dimensional system (e.g., d = 2 for a single layer). Thus, for generic d, 
V{q) = Ae^/|q|''~^, where A is a constant which depends on the dimension d, e.g., A = ^tt/eb for d = 3 (isotropic 
three-dimensional system) and A = 27r/eB for d = 2 (single layer). Here es is the dielectric constant of the ionic 
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background. The background ensures overall charge neutrality and cancels out the apparent divergence of the Coulomb 
interaction V{q) at q = 0. We point out that Eq. (|13|1 generically describes a density-density interaction mediated 
by the potential V^(q)- The particular choice of ^(q) distinguishes the short-range from the long-range case. 

After including the term Ijl^l) in Eq. ||SJ), wc follow the standard procedure and introduce a Hubbard-Stratonovich 
field phs associated to the electron density, thus decoupling the Coulomb interaction in the particle-hole channel. 
Since the short-range pairing in Eq. Q is important for small center-of-mass momentum, while the Coulomb effects 
are important for small momentum transfer, one can reasonably assume that the breakup of the actual interaction 
in this manner is physically sensible and does not lead to any "overcounting" . Note that in principle also the pairing 
term should be decoupled in the particle-hole channel, but these terms are safely negligible, at least for weak and 
intermediate U. Nevertheless, in Sec. IVll we shall discuss the subtleties of such a procedure, in the case of strong 
pairing interaction. 

As described in App. ^ by separating pHsio) = Phs{0)+Sphs{q), and by neglecting the fluctuations of the modulus 
(5|A|, the expansion of the phase-only action around the BCS saddle point leads to Scharged = Smf + Seff{9, pus), 
where now 



S, 



eff[ 



ISpHs) = X! 



2V{q) 
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(14) 



and 
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The Hubbard-Stratonovich term Spjjg{0)/V{0) does not appear in the mean-field action Smf as it is supposed 
to be canceled out, in order to preserve the charge neutrality of the system, by an equal (and opposite in sign) 
contribution coming from the interaction with the ionic background. This term survives instead in the case of short- 
range interaction. Since the self-energy Sc depends now also on the second Hubbard-Stratonovich field Spns, to obtain 
the phase-only action, one has to proceed in two steps, as is explicitly shown in App. first, the Gaussian action 
in both 9 and Spns is derived and, then, after integrating out the field Spnsi one obtains the Gaussian phase-only 
action 



(16) 



which generalizes Eq. H10() for a charged system. The coefficients £, whose definitions are reported in App^ are 
simply the coefficients A of the unscreened phase-only action ()10|l dressed by the density fluctuations within the 
random-phase approximation (RPA). This represents the only difference between the Gaussian actions Hl()|l and H16|l . 
which have therefore the same formal structure. 

The inclusion of Coulomb interaction obviously reflects on the collective mode of the system. Indeed, the expression 
of the Gaussian action H16|l in the hydrodynamic limit differs from Eq. 1)11(1 only in the coefficient of the term, 
since both A^^ and vanish in the static limit, and £}j behaves as A}j, whereas Cpp gives the RPA compressibility 
of the charged system 
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CpM, - 0) ^ ^ ^ y(q)«;o(T) ^ y(q) ' 



which vanishes for q ^ 0. As a consequence, the collective mode described by the action 
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e{q)e{-q), 



(17) 



(18) 



is the plasma mode in d dimensions. For example, in the three-dimensional case Eq. (I18|l reads 



oG 

^charged 



q 



Aire'' 



D{T) 



(19) 



and the plasma mode is characterized by a finite energy at zero momentum, Q.p — A-Ke^DiT) — ^ Ane'^ps (T) / m. 
At T = 0, where ps coincides with the particle density, D,p is exactly the plasma frequency of a three-dimensional 
charged system. 
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The vanishing of the compressibihty (|17() in the presence of density-density interactions is pecuhar of the Coulomb 
case. If one considers a short-range interaction V{q) = V, Eq. H17|l simply gives the RPA dressing of the bare 
compressibility 



1 + Vko{T)' 



C,p{q,n^ = 0)^K(T)^^^^^, (20) 



which inserted into Eq. p6(l leads to the Gaussian action 



s: 



G 
hort 



r"m "T 1 



0{q)e{-q). 



1 + Vko{T) m 

In the presence of short-range density-density interactions the collective mode is still a sound mode (known as 



Bogoljubov- Anderson modc^°), with velocity = ^/pjT^/mMT) = Csdl + Vkq{T) 



III. THE EFFECTIVE ACTION FOR THE LATTICE MODEL 

We extend the standard procedures discussed in Sec. ^ to deal with the more realistic case of a lattice model. This 
extension does not reduce to a trivial modifications of the continuum results, the most important differences being 
the appearance of vertices with more than two incoming dashed lines, and the impossibility to absorb the dressing of 
the fermionic lines in a chemical potential shift. 

We first consider the neutral case, the generalization to the charged case being straightforward. The main difference 
between the continuum and the lattice model is in the kinetic term of the Hamiltonian. Here, we rewrite the free 
action (see Eqs. Q, ((SJ) introducing a hopping t between nearest- neighboring sites on a lattice of spacing a = 1, 

S^^cro = f drlY, ctAr){dr - M)cx.(r) ~ t J2 i4Ar)c.'Ar) + h.c] + Hi{r) \ , (21) 

so that the free-electron dispersion now reads — J^j cos kj — /i, with j = xi , . . . , x^. The different form of the 
kinetic term modifies the effect of the gauge transformation on the free action Sq. As we showed in Sec. ^ after 
the transformation ljA2p the phase disappears from Si and modifies the free action 5*0 according to Eq. ©. As a 
consequence, to obtain the self-energy E within the lattice model we analyze the effect of the gauge transformation 
on the free action (|21|l . the subsequent steps being unchanged. We find two differences with respect to the continuum 
case: 

(i) The terms 'Sp and of the self-energy lO arise, after the gauge transformation, from the kinetic term V^/2m 
of the free action. Since the analogous of the velocity k/m of the continuum case is 9^k/9k in the lattice case, 
it follows that for a lattice model each insertion of a spatial derivative of in the fermionic loops is associated 
to a same-order k-derivative of the band dispersion ^k- 

(ii) On the lattice the coupling of the fermions to 9 induced by the gauge transformation of the kinetic term generates 
higher than second-order derivatives of in the self-energy E. In the language of Feynman diagrams this means 
that besides the vertices with one or two incoming 6'-lines depicted in Fig. ^ we should add all the vertices 
with n > 3 incoming 6'-lines, each associated with a n-order k-derivative of fk, see Fig. 13 Observe that no 
mixing between different components {a = x,y, ■ ■ •) of the k-derivative of ^k is produced, at least in the case of 
nearest-neighbor hopping. 

As a consequence of (i) and (ii), the self-energy E(fc) appearing in Eq. © for Seff{0) reads now 

a=xi...i:j 

0iq,) . . . 0{q^)S [J2q^^k' + k)-^ f (n) (22) 



where S{q) is a delta-function and f (n) corresponds to tq for n odd (fermionic-like contributions) and to 73 for n even 
(bosonic-like contributions). The functions sin(qi/2) in Eq. (|22|l reduce to the factors qi/2 of the continuum case at 
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FIG. 3: Vertices arising from the kinetic term of the lattice model after the gauge transformation. Each vertex with n 9-Vmes 
is associated to a same-order derivative of the band dispersion. The first diagram is present also in the continuum case (see 
diagram (b) in Fig. 0, with a factor 1/m. 



small momenta, while preserving the lattice periodicity. Since we are interested in deriving the hydrodynamic action, 
in the following they will be approximated by qi/2. 

In this section we focus on the generalization of the Gaussian action derived for the continuum model to the lattice 
case. To this purpose, it is sufficient to take into account the modifications described in the item (i). Indeed, at 
Gaussian level, only fermionic loops with two 9 lines contribute, i.e., the same diagrams depicted in Fig. |21 which 
acquire however a different meaning. In particular, the first fermionic loop of Fig. 12 which corresponds to the 
diamagnetic contribution to the current, does not correspond to the particle density divided by the mass, resulting 
instead (in the specific case of a square lattice with nearest-neighbor hopping) proportional to the average kinetic 
energy of the particles 




where we indicate with the superconducting quasi-particle energy defined in App. A, and, for the sake of clarity, 
we made the dependence on the lattice spacing a explicit. As a consequence, the current-current correlation function 
Ajj becomes now Ajj — —TSat + ^fj{l)i where it is assumed that in the definitions l|A12|l - ljA14|l each factor k^/m 
is substituted with 9^k/<9ka and the lattice dispersion appears now in the quasiparticle energy E\^. The stiffness, 
defined as the static limit of the transverse part of Ajj, is the lattice analogous of Eq. H12|l . i.e. 



D{T) = 




At Gaussian level the only difference in the phase-only action between the continuum and the lattice case is in the 
definitions H12|l or H24|l of the superfluid stiffness, the remaining steps being the same. In the following, for the sake 
of simplicity, we write down the generic Gaussian actions as 

= ^ E [^"™ + ^(^)q'] ^('Z)^(-'Z), (25) 
q 

where x indicates the static limit of the density-density correlation function xW) = ^ppil)^ ^^"^ compressibility, 
for neutral systems 

^(9) = -'^Y.^AQo^k)T3Qo{k + q)h], (26) 

k 

^ = ly [4tanh^-2/'(i?k)4 

which has to be replaced by XLR{^) ^ 1/^(q) in the presence of Coulomb interaction. As a consequence, in the 
hydrodynamic regime, the Gaussian propagator for phase fluctuations will be generically given by 

P{q) = {9{q)e{-q)) = 4 [xiq)^^ + Dq'] . (28) 

The stiffness which appears in the Gaussian propagator is the bare one, which we indicated with a subscript Dq 
in Sec. ^ However, to simplify the notation, in the following we drop the subscript, while keeping in mind that the 
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Gaussian value D is corrected by anharmonic phase fluctuation according to D ^ D + SD, where 6D will be calculated 
within the one-loop approximation. When quantum corrections are finite, i.e. SD{T = 0) ^ 0, the zero-temperature 
stiffness does not coincide with the Gaussian value, defined by the T = limit of the BCS formulas H12|l and H24|l . 
The consequences of such an occurrence will be discussed in Sec. IIVI 



IV. ONE-LOOP CORRECTIONS TO D WITHIN THE CONTINUUM MODEL 



Before deriving the anharmonic terms in the phase-only action, we discuss the phase-fluctuation correction to the 
stiffness at T = in the quantum XY model, which was briefly anticipated in Sec. ^ The quantum XY model 
includes quantum effects for the phase through a suitable modification of the Gaussian phase propagator only, while 
the interaction terms are powers of (di — dj) ^ Co|V0|, as in the expansion Consequently SD = — (_D^Q/2(i)((V0)^), 
which generalizes Eq. in d dimensions. Since ((V0)^) = q'^P{q), with the quantum propagator P{q) defined in 
Eq. (Ei), we find 



where b{x) = [e^ — 1]^"'^ is the Bose function, Cg = \J D jx is the sound velocity, g = AdC^^^ I d{d + 1), Ad is the solid 
angle in d dimensions and C ~ 1/Co is the momentum cutoff.^^ We considered the neutral case as an example, the 
generalization to the charged case being straightforward. 

To compare this result with that expected for the microscopically derived phase-only action, a preliminary distinction 
arises between continuum and lattice models. The Gaussian phase propagator is given by Eq. H25|l in both cases, 
and the temperature dependence of D{T) is due to quasiparticles only. In spite of this analogy, the stiffness D has a 
different meaning in the continuum or lattice case. In the former D = Ps/tu^ with the constraint p = ps + Pn at any 
T, whereas in the latter D is related to the mean kinetic energy, see Eqs. (fT^ and (|^ . respectively. In particular, 
at T = 0, 

continuum model ^ 0(0)^ = — , (30) 



m 



lattice model 



Z,(0,.T,0).±|:(||)(l-|)^ ,31, 



The equality ps — p is expected to hold at T = in the continuum model, since Galilean invariance ensures that at 
T = 0, where quasiparticle excitations are absent, there is no preferred reference frame in the system, and all the 
electrons contribute to the supcrfiuid component. Thus, Eq. I|30|) states that the Gaussian approximations preserves 
that Galilean invariance of the microscopic model. Furthermore, despite the fact that, beyond the Gaussian level, 
anharmonic terms induce corrections to the bare Gaussian stiffness, such corrections are expected to vanish at T = 
in the continuum model, to preserve the Galilean-invariant relation Ps{0) — p- In the lattice case, instead, the presence 
of finite corrections to the Gaussian value D{0) in Eq. (|31|) is not forbidden. This is exactly what happens in the 
quantum XY model, which is a lattice model, according to Eq. H29|l . As a consequence, the quantum XY model 
cannot be a good approximation for the phase-only model of the continuum system, since it misses the cancellation of 
the T = correction to D due to anharmonic phase fluctuations. This requirement is instead fulfilled by expanding 
the microscopically derived Seff{0) beyond the Gaussian level, and deriving one-loop corrections to Z?, as we shall 
show in the following. 

The failure of the quantum XY model for continuum systems is not a purely formal question, as it rises the issue of 
its quantitative validity both in the lattice case at T = 0, and in lattice and continuum systems in the classical regime, 
above some crossover temperature Td (which depends on the parameter of the microscopic model). This issue is of 
particular relevance when high-Tc superconductors are considered. Indeed, although in charged systems Td can be 
a quite large energy scale, of the order of the plasma energy, dissipation in d-wave superconductors strongly reduces 
Tci , as discussed in Ref. jGi] , making the classical regime experimentally accessible at much lower temperature. This 
situation then requires a critical discussion of the corresponding classical model for phase fluctuations, with respect 
to the standard assumption of considering the XY model. Here we do not discuss neither dissipative effects, nor the 
precise deflnition of T^i, but we directly identify the classical regime by considering only the Vim — contribution in 
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the Matsubara sums which define the phase-fluctuation corrections 5D and SDxy- In the case of the XY model, by 
considering Eq. H29|) . wc find that in the classical regime 

SDxY « -2-^, (32) 

i.e., SDxy{T) oc T, as already observed in Sec. |l| We shall find below that also 6D is linear in T, but with a much 
smaller coeSicient in the weak-to-intermediate-coupling regime (see also Ref. 0). 



A. Classification of the diagrams 



We start our analysis with the continuum neutral system, while the long-range Coulomb forces will be included 
later. To determine the one-loop correction SD to the stiffness induced by anharmonic terms in Seff{0), we derive the 
contributions up to ^4 in Eq. @. In general, if we separate the effective action in the Gaussian (5 ) and anharmonic 
{Sank) part, and make a Taylor expansion of the latter, by means of the Wick theorem we can evaluate the two-field 
propagator as the resummation of all the connected diagrams obtained by fixing two 9 legs and contracting the others 
with the Gaussian propagator P{q), 



{9{q)0{-q)) = 



/g(g)0(-g)e-^''(^)e-^°"'-Wp0 

^-^ n] Z ^-^ TV. 



where {•■■)con indicates connected diagrams only. All one-loop corrections to P{q), i.e., diagrams with only one 
integration on bosonic momenta, are obtained by approximating X)n(~l)"'/''^'['5'anft]" — S'4 -I- [83]'^ /2. Moreover, since 
we are specifically interested in the correction to D, which is the coefficient of the gradient term in the effective 
Gaussian action, we select the diagrams which have at least two \79 lines, as depicted in Figs. 0]and[Sl 





3a 36 

FIG. 4: Vertices of S3 which contribute to SD. We use the same notation of Fig. Q 
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FIG. 5: Vertices of S4 which contribute to SD. According to our classification, 4a is a fermionic vertex, 46, 4c, 4d are mixed 
vertices, and 4e, 4/ are bosonic vertices. 



To obtain one-loop corrections to D we consider all the possible contractions of two 9 lines in one of the fourth-order 
vertices in Fig. [3 or in the combination of two third-order vertices in Fig. 0] The above classification in bosonic, 
fermionic and mixed contributions turns out to be particularly convenient while evaluating all the resulting one-loop 
corrections to D. We can indeed distinguish: 
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(Z31) (B2) (S3) 



FIG. 6: Bosonic corrections to D. Bl arises from the diagram 4e of Fig. |S] B2 and S3 from two 36 vertices of Fig. 3] Notice 
that the dashed-dotted line of S3 indicates the average {Qrn<i0{<l)0{~q)) ■ 

• Fermionic corrections, which are given by the self-energy- and vertex-like diagrams obtained from the 4a loop 
in Fig. El 

• Bosonic corrections, which arise from two loops of the 35 type and from the loop 4e of Fig. El see Fig. El 

• Mixed corrections, which are the remaining combinations of third- and fourth-order diagrams. 

The distinction between fermionic, bosonic, and mixed contributions to SD is motivated a posteriori by the fact 
that all the contributions to SD coming from the fermionic and mixed diagrams cancel out each other at T = 0, 
whatever is the symmetry of the gap and both in continuum and lattice models. This result is connected to the 
fact that the temperature dependence of the fermionic and mixed bubbles is controlled by the thermal excitation 
of the quasiparticles, which are suppressed at low temperature. However, while for s-wave superconductors this 
suppression is exponential in temperature due to the presence of a finite gap in the excitation spectrum, for a d-wave 
superconductor the presence of nodal quasiparticles obliges a careful evaluation of the fermionic contribution to SD. 
For example, one can see that the fermionic contributions generated by the vertex (4a) in Fig. El can be written as 

'^^e = E ^('z)q'Q(q,a„), (33) 

where vi^ — 9^k/9k and the hydrodynamic hmit of the Q bubble is 

Q{ci^Q,^^ = Q) = ^Y.'ir{E^), (34) 

k 

where f"'(x) is the third-order derivative of the Fermi function. In the s-wave case the integral in Eq. (|34|l vanishes 
exponentially at T = 0, while in the d-wave case, by evaluating the contribution of nodal quasiparticles to Eq. I|34|l . 
which is the relevant contribution at low temperature, one finds that it diverges as l/T. This can be easily seen by 
considering that the density of states of nodal quasiparticles is linear in energy N{E) oc so that Q{q = 0) oc 
/ dEN{E)f"'{E) l/T as T ^ 0. It should be stressed that the divergence in the d-wave case of all fermionic 
loops with > 4 Green functions, when evaluated at zero external momenta, is a characteristic feature of systems 
having Dirac quasiparticles. Nevertheless the one- loop expansion of D is well defined: indeed by retaining all the 
q, Qrn dependence on the fermionic bubble Q one has no longer to deal with any divergence of fermionic corrections 
to D, which result even vanishing in the zero-temperature limit. 

Special attention must be devoted also to the bosonic (dl and d2) and fermionic (d3) diagrams depicted in Fig. 
[3 Indeed, these diagrams appear as corrections to the bare Green function in the diagrams of Fig. (21 which define 
the the mean- field current-current correlation function Ajj in Eq. (|A13|) . In particular, dl and d2 in Fig. [71 are 
phase-fluctuation corrections to the diamagnetic term — pmf/w, and d3 in Fig. [71 is a phase- fluctuation correction to 
the bubble Ajj, which gives the quasiparticle depletion of Ps(T) at finite temperature. One can easily see that these 
terms are due to the shift of p with respect to the mean- field value pmf, as induced by anharmonic phase fluctuations. 
On the other hand, working at fixed density, one assumes pmf = P and takes into account these contributions by 
shifting the chemical potential I-J, I-J, + Sp with respect to the mean-field value, so that 



(35) 
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where we use the definition x = dp/dfx for the compressibility, and the free-energy contribution AF due to phase 
fluctuations is evaluated within the Gaussian approximation, 

AF= -In J P^e-^'^W = 1^ ^ H^fl^ + Dq^]. (36) 

While evaluating one- loop corrections to D, one takes into account the chemical potential shift H35(l by adding the 
contribution 

SD, - ^Sf,, (37) 

which cancels out exactly the contribution of the diagrams in Fig. {Tj). This can be easily understood, e.g., at T = 0. 
Because at T = all the mixed diagrams vanish, only dl and d2 contribute to SD. Let us consider the diagram d2. 
Its contribution to 6D is 



1 T 



q (. k ) q 



Indeed, the fermionic loop of d2, given by the trace in the previous equation, corresponds to the hydrodynamic limit 
of the same density-density bubble x{l) which appears as the coefficient of the term in the Gaussian propagator 
P{q) (EHl, see Eqs.^-lISZIl-^^ 

On the other hand, while evaluating d/S.F/d^ according to Eq. (|36|l . one finds that the Dc^ term of Eq. 13()|) 
contributes with 

dAF , ir>r^ 1 J_Lsr 2pr 



d^i x^l, + Dci^ d^i 8m r2 
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where we used Eq. and D = ps/m = [p — p„(T)]/m so that, a.i T = D — p/m and dD/dp — xl"^- As a 

consequence, by means of Eqs. (|35() and H37|l we find that 
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which is exactly —5Dd2- The same holds for SDdi, SDds and the contribution arising by differentiating x^m ^-^^ 
Pn{T) (at finite temperature) in Eq. 1)36(1 . leading to a cancellation of (51?^ with the diagrams of Fig. [3 






(dl) (d2) (d3) 



FIG. 7: Bosonic (dl and d2) and mixed (d3) one-loop corrections to D arising from phase-fluctuation correction of the electron 
density beyond Gaussian level. Working at fixed density, these terms are canceled out by the SDfj, contribution due to the 
chemical-potential shift. 



B. One-loop corrections at T = 

The cancellation of fermionic and mixed contributions as T ^ can be demonstrated with lengthly but straight- 
forward calculations, so we will not discuss it here, where instead we analyze the behavior of the remaining bosonic 
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contributions of Fig. |B1 and in particular the role of dynamics in restoring Galilean invariance in the continuum 
model. We note in passing that the cancellation of the fermionic corrections has been already observed in Ref. [l^ . 
where however the authors did not include both mixed terms and corrections to 5D coming from the third-order loops 
in S^{9). While neglecting mixed term at T = is not problematic, since they cancel out at T = 0, neglecting the 
dynamic third-order vertex 36 of Fig. ^ leads to neglecting the bosonic diagrams B2 and Bi of Fig. (BJ whose role at 
T = is crucial. 

Roughly speaking, the diagrams depicted in Fig. |^ arising from the bosonic self-energy Y,b only, correspond to 

2 

B\ ^ -^(vefuvef), B2, B'd ^ ^(vefuvefidref), 

where the angular brackets indicate the average with respect to the Gaussian action. As a consequence, the first 
diagram Bl, arising from a vertex ~ in the effective action S^fj, is the analogous of the term J2a=ii^ctO)^ in 

the quantum XY model Indeed, the term Bl generates the correction 

^^-^-4il2 5:^^^('^M^)' (38) 

9 

where xil) is given by Eq. H26|l . By evaluating Eq. H38() in the hydrodynamic limit, in analogy with Eq. H29(l . we 
find that 

= ^^^^^ - E 7^[^(^^^'^') - '^-'^^^'^'^l ^ - Tfe 2^' ^''^ 

which is again a finite correction, as in the quantum XY model, see Eq. H29|> . By comparing Eq. H29|l and 139|l . one 
finds that the strength of the term (V6')'* in the continuum model is controlled by x/m^, and not £,qD, as in the XY 
model. Nevertheless, at T = 0, the result (|39() cannot be conclusive, since it would lead to ps ^ p at zero temperature. 
Indeed, by adding the contribution of the remaining diagrams B2 and Bi 

q 



we find that the bosonic one-loop corrections to D reduce to 

1 T 



^^ = - 7) E q'^(9)x(9) [2 - nlP{q)x{q)] . (40) 
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In writing Eq. H4()(l we are relying on the fact that both the third- and fourth-order vertices needed to calculate 
SD are expressed in terms of the fermionic bubble x{q)i Eq. H26() . Indeed, both the static — ((V6')^/8m)f3 and 
the dynamic E]j = {idrO/'i)t^ part of the "bosonic" self-energy are tied to the Pauli matrix fa, leading to the 
same fermionic bubble (|26|l as a coefficient of the phase-only action. As observed in Refs. J.9j , this is the relevant 
consequence of the Galilean-invariant form of the bosonic ra term appearing in Q. Evaluating Eq. H4()|) in the 
hydrodynamic limit, we find 



6D 



2 



2m? dn ^ x^l-, + 



-nt + uj. 



1 7" 2 (»^-m + ^q)^ _ 1 r 2 ^ _ ^ L "S^ 2,// ^ z.-, ^ 

^ [-^O. , ^2 I ,,,212 - Om'i rlO 2^ ^ [_,0.__ ^ „, J2 ^ 0^2^0 2^^^ "I'^qJ' \'^^) 



2m? dVl \-{i^mY + i^lY 2m2 d^ ^ ^ l-i^m + <^„Y ^rn? dQ 

q,0™ L V "1/ qj qf^^ L m qj ^ 



where uiq = y^D/xIql is the sound mode, and b'{x) — — /3e'^^/(e'^^ — 1)'^. As a consequence of the addition of the 
B2 and B3 contributions, the pole at ir^m — tjq < 0, responsible for a finite contribution at T = in Eq. (jlfflfl . is 
canceled in favor of a double pole at iflm — i^q> 0, leading to Eq. H41I) . where the standard Bogoljubov reduction of 
Ps/m in a superfluid bosonic system is recognizedi^i. Thus, for T ^ the integral in Eq. H41|l vanishes, i.e., by fully 
including the dynamic structure of the interaction for the phase, we obtain that in a Galilean-invariant system ps = p 
at T = 0. Notice that the result (|41|) arises from an exact compensation between static and dynamic bosonic vertices 
for phase fluctuations, both at Gaussian level, in determining the structure of the propagator P{q), and beyond the 
Gaussian level, in weighting the relative contributions of static and dynamic parts. 
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C. One-loop corrections in the classical limit 



According to Eq. H39|l . a further result of the microscopic derivation of the effective action is a reduction of the 
strength of the static interaction term (V^)* with respect to the XY model, see Eq. (|29|l . This reduction is relevant in 
the cZassica/ regime, where only the correction SD^l to D coming from the static interaction term Bl oc (x/?ti^)(V0)'* 
survives in Eq. (|^ . leading to 

SD^- ^ ^ 



which is qualitatively similar to the result of the classical XY model, Eq. (|32l) . To make a quantitative comparison, 
one can estimate the compressibility x in the weak- or intermediate-coupling regime as the value of the density of 
states at the Fermi level,^^ which, in the continuum model, is 

/dk 
■^S{eF - £k) = 2AdVi4-\ (42) 

where kp is the Fermi wave- vector. Nevertheless, also ps, i.e., p (if we neglect quasiparticle effects^l), can be defined 
through kp^ at least in the weak-to-intermediate-coupling regime (more precisely as long as Tc < Ep, i.e., U < t) 

It then follows that the ratio between the coefficient of the static interaction term (V6')^ in the continuum or in the 
XY model is 



X X 



m^Dfo mp,^2 (^^^^)2 ^ 
which means that, apart from a factor d/2, we have 



(43) 



Classical Regime => ^ — ^ . (44) 

SDxY {kp^of 

As a consequence, in the classical limit, the one-loop correction to D derived within the effective action is smaller 
than the corresponding first-order correction in the XY model as far as kp^ < ^q, i-e., in the weak- and intermediate- 
coupling regime for the pairing interaction. 

Finally, let us consider the effect of the Coulomb interaction between the electrons. As we showed in Sec. IllBI the 
presence of a density-density interaction in the microscopic Hamiltonian implies the dressing of the coefficients of the 
Gaussian action by the RPA series of the Coulomb potential V{q). As a consequence, in the hydrodynamic limit, the 
bare compressibility x is replaced by the dressed bubble xlr{(1) ~ l/^(q)j and the sound mode Wq of Eq. H41() is 
converted into the plasma mode oj^ of the d-dimensional system (see Eq. (119(1 for the three-dimensional case) . In the 
same way, when deriving the anharmonic terms in Seff-, we must now include the RPA density fiuctuations in all the 
vertices for 5*3 and S4. The one-loop corrections to D are formally identical to Eqs. l|ini) - (|mi . with x^ Xlr- Thus 
at T = we recover again the cancellations (|41|l of the bosonic diagrams, with ujq uj^. At the same time, since 
Xlr{<1) vanishes as q ^ 0, the classical (fi„ — 0) term in Eq. (gUJ, i.e., SD^i given by Eq. ((S^ . reads 

Z_1V M - ^ 

m^dDn^^^"^' m^dDXe^ 2d - 1 ' 
q 

whereas the quantum XY model leads to the same result ()32(l of the neutral case. Thus we can roughly estimate 

Classical Regime LR - • ^ £f f 



SDxY £c ykp^Q 

where ep = kp/2'm is the Fermi energy and ec = Xe^kp is a characteristic Coulomb energy scale. As a consequence, 
while within the XY model the Coulomb interaction modifies only the low-temperature behavior of p^, within the 
continuum model it affects also the high-temperature classical regime. 
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V. ONE-LOOP CORRECTIONS TO D WITHIN A LATTICE MODEL 

To extend the previous results to a lattice model, we need to take into account the discussion of Sec. IIIII in which 
we derived the phase-only action starting from a nearest-neighbor tight-binding model for electrons in a lattice. 

At Gaussian level, the item (i) of Sec. IIIII allows us to define the superfluid stiffness of a lattice model according to 
Eq. In particular, a factor ^Aa, with 

A 

appears in the fermionic bubble which defines the diamagnetic term, i.e., the coefficient of the (V0)^ term in the 
Gaussian action for phase fluctuations. As a consequence, D is controlled at T = by the mean kinetic energy T, 
Eq. (|23|l . instead of p/m of the continuum case. In deriving anharmonic terms in Seff{0), we must take into account 
item (i) in the proper deflnition of the fermionic bubbles of Figs. 0] and El which give the vertices for the interacting 
phase-only model. Moreover, we need to include also the fourth-order vertex depicted in Fig. |S1 which is peculiar of 
the lattice case, and arises from the n = A term of Eq. (|22() (see also Fig. 




FIG. 8: Fourth-order vertex of S4 peculiar of the lattice case. In evaluating the fermionic loop one includes a 9'*^k/9k* 
derivative of the band dispersion. 

At this point the calculation of one-loop correction to D follows the same steps described in Sec. IIVI The cancellation 
of fermionic and mixed one- loop corrections to I? at T = still holds, thus we can again focus on bosonic corrections 
only. As far as the corrections in Fig. dare concerned, they are no more canceled out by the contribution SD^^ due 
to the chemical-potential shift beyond the Gaussian level. In particular, while discussing in detail in Sec. IIVI the case 
of the term d2, we stressed that in the continuum model such a cancellation holds because the fermionic bubble of 
diagram d2 is the same bubble which one obtains by differentiating the diamagnetic part of D with respect to fi, so 
that dD/dfj. = x/"i- Indeed, since in the continuum case the factor A^ = 1/m associated to the (V0)^ insertion is 
a constant, one obtains the same two-line fermionic loop both by expanding Seff up to 5*4, leading to d2, and by 
differentiating the single-line diagram which defines p. However, in the lattice case the factor Aq, is not a constant. 
As a consequence the diagram obtained by differentiating the diamagnetic term T = -D(O) with respect to /i has one 
insertion of A^, while the the diagram c?2, obtained by expanding Seff according to item (i), has a A^ factor for each 
6*- vertex (diagram C2 in Fig. I^J. The same holds for the dl diagram. 

As a consequence, in the lattice case, the contributions to D come from the bosonic diagrams depicted in Fig. |51 
and from the correction SD^^ induced by the chemical-potential shift, which do not cancel anymore the diagrams Ci 
and £2 in Fig- HI While evaluating 6Dc2 we must consider that we have 2{d — 1) bubbles with A^ 7^ A^^j in the 
expansion of Sank in d dimensions, so that 

5Dci 
5Dc2 
SDc3 

6Dc5 + SDc6 



1 T 



hl^Yl ^('?)q' [XEE + {d- 1)xee] 



1 T 
8dn 

1 T 

'sdn 



1 

Y^^PiqWxEEiq), 
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FIG. 9: One-loop correction to D within the lattice model. As usual, we indicate with the dashed line the gradient of the phase 
and with the dotted line the time-derivative of the phase. Moreover, we explicitly indicate the insertions of band-dispersion 
derivatives in the fermionic loops. 



where r]pE{q) is the fermionic bubble with three Qq lines and one insertion of A^, and the Xabil) bubbles correspond 
to the insertion of one (xpfi) or two equal/different (xee/xee) factors A„ in the two-line bubble. As usual, when the 
dependence on q is not explicitly indicated, we are considering the q — > 0, 0^ = (static) limit. Notice that, since 
(9^^k/9k^) — — a^(9^^k/9k^), the T which appears in SDc^ is proportional to the mean kinetic energy defined by 
Eq. (|23|l . T — a^T = a^Ekin (the dependence on the lattice parameter a is made explicit, for the sake of clarity in 
the forthcoming discussion). 

According to Eq. H35|l . the chemical-potential shift is 

1 d^F 1 ^. ,.^2 / ^ 21 

s^^ = -^ = ^^T.Pi<i)i^lM<i) + ^\,Eh 

where ri{q) is the three- line bubble without factor A^, x is the compressibility, which carries no A^ insertion, and we 
used XpE = dD/d^i. By means of Eq. (|37|), we obtain 



5D, = ^6f^ = ^^-lj2piq)[nlrj(q)+ci\,E]. 



(50) 



Finally, the overall one- loop correction to D coming from Eqs. H45|) - (|49|l and H50|) is 

1 T 



SD= - 
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"^XEEiq) + XEE + [d- l)XEE + T - d 



xU^qWP{q) 2d^pE{q) + d^viq) 



(51) 
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To make a comparison with the continuum case, Eq. H41|l . we evaluate the above equation in the hydrodynamic 
hmit. This corresponds to calculate the Xab bubbles at zero incoming momentum and by carefully evaluating the q — > 
limit for the terms which involve r]pE{q) and 77(g). Indeed, by simply calculating X]n,„ Pil)^mV{Q) ~^ vJ2n,„ Pil)^m 
one finds that P{q)fl^ —> 00. Instead, by retaining the flm dependence in the fermionic bubble one ensures the 

convergence of the Matsubara sum, calculating the q ^ limit afterward. Thus at T = limit we obtain that the 
overall correction to D is finite in the lattice case. 
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according to the following limiting values 
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The r = value of x is given by Eq. H27() . The others bubbles are defined as 
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(53) 
(54) 
(55) 

(56) 
(57) 



Thus in the lattice case the T = mean-field value of the superfluid density is corrected by one-loop phase 
fluctuations according to Eq. (|52|l . It is interesting to observe that the continuum case is recovered from Eq. (|52|) by 
performing the limit a — s- 0, i — > 00, keeping 2a^t = 1/m finite. Indeed, since = — 2iek — /i, with = cos(afca), 
Aq — {d'^£,k/dkl) = 2ta'^ek ~ 1/m in this limit, and according to Eqs. (|^ - lf57|l . one finds xee,Xee x/^'^: 
XpE — > x/iTT-j VpE ^ ry/m, while T = a^D — > 0. Thus in this limit, Eq. (|52|l correctly reduces to the contimmm 
(Galilean- invariant) result SD = of Eq. (|41|l . 

We now compare SD and SDxy in the lattice case. Differently from the continuum case, this comparison is now 
meaningful both in the classical limit and at T = 0, since now we find finite corrections to D. By comparing Eq. 
(|52|l and Eq. (|29|) at T = 0, we see that roughly speaking the strength ^ ^qD of the interaction in the quantum XY 
model is given by the quantity in square brackets in the first line of Eq. H52|l . i.e., it is of order 
Following the same arguments which led us to Eq. 14311 . we can estimate 



Xee,Xee,X^peIx- 



XEE, XEE: 



X^pE 



|2 ' 



where XeEtXee,X^peIx have the same role as x/to^ in Eq. H43I) . and kp ~ 1/a in the lattice case. However, in the 
lattice case the dynamic contribution controlled by XpeIx^ does not cancel with the static contribution, and both 
contribute to SD at T = 0. Analogously, one finds that T = a^D « D/k"^ = i'^D / ikpiaY ■ Moreover, 
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where N{eF) is evaluated in Eq. H42() . A is the maximum gap value, and we used the fact that at weak coupling 
kp£,Q ~ kpVp/A — fi/ A ^ Im^ Taking into account the prefactor g/Cx in the second line of Eq. (|52|l . we have 
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Cx 



XpE ~ 

VpE V 

X 



6DxYj^^ (58) 



where f/pp ~ rj/m. Thus, the contributions coming from the three- line fermionic bubbles are subleading, and at T = 
Eq. leads to an estimate 



SD 



SDxY {kp^o 



.2 ■ 



(59) 



of the same order of Eq. I|44(l . within a numerical factor. The presence of Coulomb forces does not change qualitatively 
this conclusion, and introduces only minor quantitative corrections. This is due to the fact that the RPA expressions 
for xee, Xee have a finite limit for q — s- 0, contrary to the continuum case, 

LR_^ xUVj^) ^, / xIe\ xIeH'-' 
Xee - Xee ~ , , ... x - Xee H r-j-- (oUj 

When (5D is evaluated with the long-range phase propagator /*((?), one sees that the leading terms in powers of 
l/(CofcF) come from the finite q = limit of xee-, Xee in Eq. H6UI) . and from the T term, all of order l/(^ofcF)^ with 
respect to the coefficient of the XY model, whereas 

LR ^ XpE ^ 1 

i + y(q)xpF 'i^o nq)' 

so that dynamical terms are subleading at T = in the long-range case. This is clearly true as far as the constant 
in the parenthesis of Eq. (|60|l is a finite quantity. Indeed, when it vanishes by approaching the continuum limit, all 
(static and dynamics) terms are of the same order, canceling each other as discussed above. 

Finally, the same result H59|l holds in the classical limit, both for the neutral and for the charged system. As we 
discussed above, at high temperature only the static correction to D contributes, with a coefficient controlled in the 
neutral case by XEE^XEE^^' ^ or by their q = finite limit in the charged case, which play the same role (and have 
the same estimate) as x/m^ in the continuum case. 

As far as the effect of bosonic degrees of freedom is concerned, our result H59|l suggests that the XY model can 
eventually lead to an overestimate of phase-fluctuation effects with respect to the result obtained by means of the 
microscopically derived effective action. Indeed, according to Eq. (|59|l the depletion of D due to phase fluctuations 
can be quite small at weak and intermediate coupling (particularly in the BCS limit) both in the quantum and in the 
classical regime. 

On the other hand, as far as the temperature dependence of 8D for temperatures < T < Tc! is considered, one 
should take into account also the effects of fermionic and mixed contributions. As we discussed above, quasiparticles 
contribute to the depletion of D at finite temperature already at Gaussian level, with a temperature dependence 
determined by the symmetry of the order parameter, which is linear in T in d-wave superconductors. One would 
expect that at T > the one-loop corrections to D coming from fermionic and mixed diagrams correct quantitatively 
but not qualitatively the D{T) dependence due to quasiparticles at Gaussian level. However, a detailed analysis of the 
temperature dependence of D(T) as arising from bosonic, mixed and fermionic contribution at low temperature is a 
much harder task. From one side, one has to face the difficulty of evaluating the q = limit of these expressions, as we 
discussed in Sec. IV A. On the other side, at finite temperature, non-local terms appear when deriving the phase-only 
action, making the hydrodynamic expansion not well-defined. These are the so-called Landau singularities, arising 
from the non analytic behavior at small q of the AjjpH In any case the analysis of these issues requires separate 
investigation and is not addressed in the present work. 



VI. THE STRONG-COUPLING REGIME 



The results for the lattice case was derived by assuming l/{kp^o) as a small expansion parameter. In particular, 
this assumption allowed us to estimate the fermionic bubble according to Eq. 1)4211 . i.e., x — ^{^f), and to neglect 
the contribution to SD of the fjab bubbles, see Eq. (|58|l . One should expect that such approximations are reliable 
in the weak- or intermediate-coupling regime, where the correlation length is larger than the lattice parameter 
a, which sets the scale of the Fermi wave- vector kp ~ 1/a. However, in the strong-coupling regime for the pairing 
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interaction several modifications should be included. The fermionic bubbles can no longer be estimated as the density 
of states at the Fermi level; the correlation length is generically expected to be of the same order of magnitude as 
the lattice parameter; one should incl ude both the RPA fluctuations, also induced by U , in the particle-hole channel, 
and the fluctuation of |A| ~ {U /2)yJ p{2 — p), which fluctuates because p fluctuates.'' This last issue can be better 
understood by analogy with the superfluid bosonic systems, where the order parameter is directly the square-root of 
the density, and one expects that, in the strong-coupling regime, the fermionic system can be mapped onto an effective 
bosonic system, where Cooper pairs act as bosons with a weak residual repulsion. The above points are intimately 
connected, as we already discussed in Ref. ITsl There, the behavior of the correlation length, properly defined as the 
characteristic length scale for the fluctuations of the modulus of the superconducting order parameter, was analyzed 
in the strong-coupling regime. The main finding is that attains a value of the order of the lattice parameter a 
at finite densities, and diverges as ^/^/p in the low-density regime p ^ 0, recovering the well-known behavior of a 
bosonic systemi^ 

In the following we consider the s-wave case, which allows us for a transparent analytical treatment. We first 
address the neutral case, i.e. the model Q-Q with w(k) = 1 (negative-f/ Hubbard model), and then we extend our 
results to the charged case. In the presence of d-wave pairing the analysis becomes more difficult, even though the 
physics is not expected to be different from the s-wave case, as we briefly discuss at the end of this section. 

As we said at the beginning of Sec. Ill Al the Hubbard- Stratonovich decoupling is a useful tool which allows us 
to substitute the microscopic model with an effective action written in terms of the relevant collective degrees of 
freedom. However, the choice of the relevant variables depends on the physical properties of interest and on the 
energetic scales involved in the problem. For example, the negative-?/ Hubbard model can be decoupled both in the 
particle-particle channel, as we did, and in the particle-hole channel. The latter decoupling would lead to an effective 
action depending on the particle density as well, Seff — Seff{0, ^|A|, 6p)^ The resulting action is formally the same 
that we obtained by decoupling the additional interaction Hamiltonian (|13|l in the particle-hole channel, but with 
V{q) replaced by the same —U/2 which appears in the Hubbard term. The relevance of density fluctuations within 
the Hubbard model at strong coupling is addressed in Refs. In deriving Seff up to Gaussian level a term 

proportional to the product of amplitude and density fluctuation, ^ (5| A|(5p, appears. The coefficient of such a term in 
the Gaussian action is vanishingly small at weak coupling. This means that density and modulus fiuctuations decouple 
at weak and intermediate coupling, and one can safely neglect modulus fluctuations without loosing information on 
the contribution of the density fluctuations. As a consequence, one can derive the phase-only action according to 
the procedure described in the previous sections. However, at strong coupling the situation is more involved: as it is 
shown in Refs. j by increasing U the density- modulus coupling become sizable, and one finds that the density 

of particles and the modulus of the order parameter experience the same fluctuations. Since the phase mode has 
the same behavior of the density mode, in the strong-coupling regime we should also take into account amplitude 
fluctuations to recover a consistent description of the phase fluctuations. 

In Sec lIIBI we found that in deriving the effective action for the phase of the order parameter the inclusion of density 
fluctuations reflects in the replacement of the "bare" coefficients A of the Gaussian action with their counterparts 
£ evaluated within RPA resummation in the particle-hole channel. Beyond Gaussian level one expects the RPA 
dressing of the various fermionic loops which appear as coefficients of the effective action. As far as the 6D correction 
in Eq. (|52|l is concerned, this leads, e.g., to substitute the Xab fermionic bubbles with the RPA resummation with the 
potential —U/2 of the corresponding density or current irreducible correlation functions X^abil)- then left with 

the problem of including modulus fluctuations. In principle, to take them into account one should derive the Gaussian 
effective action Sg{0, S\A\,Sp) by including all the fields and then integrate out modulus fluctuations as well, see, e.g., 
Refs. 0^3- However, the contribution of amplitude fluctuations to the static long- wavelength limit of the fermionic 
bubbles X^Ibi^ ~^ 0,t^n = 0) can be easily derived in a different way. Let us consider, e.g., the density-density 
correlation function x^^^ ■ Since x"''' = dp/ dp, (see also Sec. IIV|I . we can deduce it from the mean-field expression 
(jA15|) of the density of particles, by taking into account the dependence on p of the gap amplitude |A|(0) = Aq 
itself, which appears in the definition of i?k = Vsk+^o (^^^ ^^^'^ ^<is. (|A10|l V At weak coupling dA^/dp can be 
neglected, and one finds the definition (|27|l of x^^^ ■ At strong coupling, where modulus and density fluctuations are 
proportional, by working at fixed number of particles the chemical-potential variations needed to preserve the particle 
number reflect into amplitude variations, and dAo/dp gives a significant contribution to the irreducible bubbles. We 
then have (at T = 0): 

„r ^ dpMF ^ 1 ^ , 1 ik . 9Ao 

^ dp n^Ei n^Ei^'dp' 

where 9Ao / dp can be obtained from the self-consistency saddle-point equation 
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by differentiating botfi sides witfi respect to /z, 
so tfiat 

, 1 [Ek(ek/i?^)] 



irr 
A 



In the limit of strong pairing interaction, we can evaluate the previous bubble for t/U <^ 1, by taking into account 
the strong-coupling solution of the coupled saddle-point equations for the amplitude Aq Eq. (|61|l and for the particle 
number Eq. IjAlSjl . 

M = -|<5(l + 2da2) + o(^^ 
where we introduced the short-hand notation (5 = 1 — p and the small parameter a = 2t/U . We then obtain 

X'" = I (1 - ida") . 

Thus, including also the RPA resummation of x*'"'' in the particlc-holc channel, we finally get 

^ " T~~(C772)^ " 2dM' ^^^^ 

The expression ()fil-t(l for the density-density electronic correlation function at g = is the same which one would 
obtain by explicitly integrating out both the 5p and the i5|A| field in the action S{9,5\/!^\,8p). The latter procedure 
was indeed adopted in Ref. [13 while deriving the propagator for the Hubbard-Stratonovich field associated to the 
density. The result H()3I) shows a large increase of x at large U: this could be expected, since the electrons are strongly 
paired, naturally increasing the compressibility of the system. 

As far as the other fermionic bubbles are concerned, we follow the same steps. In analogy with the previous 
discussion 

_ ^ lao 



where the stiffness is defined in Eq. (|24() . and 

Dee = 



a,t T — 0. By taking the derivative with respect to the hopping t one includes also dAo/dt, as obtained again from 
Eq. (|61l) . In the strong-coupling limit 

x;'S - -4^2,5, (64) 

xk^E; = 2tail-S'), (65) 

xVe = 0{ta'). (66) 

By introducing the dressed potential U — (i7/2)/[l — {U/2)x^^^] ~ U/Ma^, in analogy with Eq. (jHSl, we find 

xp£ = x;i[i+x;7^] = -^' (67) 
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8ta 
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Xee = xVe + {xTeTU = 2to(l - + _52, (68) 



XEE = Tee + {.xTeTV = —b\ (69) 
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The inclusion of the modulus fluctuations always adds a second term to the weak-coupling definition of the irreducible 
bubbles, see, e.g., Eq. (|62|l . This additional term is never of lower order in a with respect to the term obtained with 
the weak-coupling definitions 1)53(1 - 1155(1 . Since one expects that the same holds for the rj bubbles, let us first estimate 
the RPA resummation of the weak-coupling definitions ((56(1 and ((57(1 . which do not include modulus-fluctuation 
contribution. The resulting corrections to D, arising from the second line of Eq. ((52(1 . are subleading with respect to 
the contribution coming from the Xab bubbles, first line of Eq. ((52(1 . Since modulus-fluctuation contributions do not 
change this result, we can neglect the correction SD coming from the r/ bubbles. 
By using the expressions (|^ - (|^ . the strong-coupling expression 

D = a^~'^2ta{l - 6^), 

and the fact that T = a^D, we find that Eq. ((52() . at large U/t, reduces to 

6D = ^8a4-'*to(l - S^) = ^a^D, (70) 

where we reintroduced explicitly the dependence on the lattice parameter a, to allow one for a direct comparison with 
the result of the XY model, Eq. (|^ . We then obtain that 

5D^5DxyI^£^ ■ (71) 

By taking into account the previous results for ^or^ we find that, since at finite densities ^ ctj the XY model 
appears as the proper "strong-coupling" effective quantum model for phase fluctuations in a lattice. A noticeable 
exception is found in the extreme low-density regime, where diverges as l/^fp. Moreover, since g ^ Eq. 
((70(1 shows that in the low-density limit the relative superfluid-stiffness correction 5D/D vanishes as p'^^'^. Observe 
that instead the XY model would lead to a relative correction 5Dxy /D ^ p'^/^-i which does not vanish at p = in 
d = 2. 

In the classical regime, as we already discussed in the weak-coupling case, only static anharmonic terms survive, 
leading to 

SD 
SDxY 

Since the quantity in parenthesis behaves as the microscopically derived (^Q/a)'^ at all densities, up to a numerical 
constant of order one, we find that the in the classical regime the XY model turns out to be appropriate, regardless of 
the density. Observe that the numerical difference comes out from the fact that the XY model attributes a somewhat 
arbitrary coefficient to the anharmonic term (V0)^. If one assumed in Eq. (^ higher order harmonics as 

^ 5](l~cos^,,)", (n>2), 
<hj> 

which do not modify the Gaussian term, by adjusting the D„ coefficients one could reproduce the correct (V6')2" 
term as it arises in the microscopic expansion. In such a case the classical limit of this extended XY model and of 
the microscopic model would lead to exactly the same result. As it has been noted recently in Refii&, the mapping of 
the classical phase-only action of a lattice systems on an extended XY model can also lead in two dimensions to an 
enhanced fluctiation region near the Kostelitz-Thouless transition. 

Let us consider now the case when also the Coulomb potential is present. The effect of modulus fluctuations on the 
irreducible fermionic bubbles is the same as already discussed above , so Eqs. ((62(1 and ((64(1 - 1(66(1 are still valid. As far 
as density fluctuations are concerned, we just have to take into account that the RPA resummation of the irreducible 
bubbles in the particle-hole channel must be performed with the full potential V{q) — U/2. As a consequence, even 
though we are considering a strong-coupling (i.e. large U) limit, as q goes to zero the Coulomb interaction V{q) is 
always predominant. Thus in the hydrodynamic limit the contributions to the correction SD only come from the static 
interaction vertices, as already discussed in the weak-coupling limit, since the XpE bubble vanishes for small q, while 
the current-current correlation functions attain a finite value. More explicitly, by performing the RPA resummations 
(EH-(|Snj with U replaced by the fuU dressed potential U -> y(q) = -[V"(q) + U/2]/[l + (F(q) - C//2)x*''1, it can be 
easily seen that 
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/ irr\2 

XEE ^ x'£'^E- -^-2H^-s') + o{ah), 

A 

( irr\2 
~irr \^pE) 3,n 

XEE Xee - = 0{a t). 

As a consequence, in Eq. H51|) for SD a coeiBcient {^Xee +'^)/^ — survives. This coefficient plays the same 

role as in the XF- model, see Eq. (|29|) (which however should be evaluated with the long-range propagator 

P{q))- Thus, Eq. (|71|l still holds, both at zero temperature and in the classical regime, where by definition only static 
anharmonic terms survive. Observe that this result is again not trivial, because only properly including modulus 
fluctuations one can find the equivalence between x^^,T and a^D in the strong-coupling regime for the pairing 
potential U. 

Finally, we briefly comment on the extension of the previous results to the case of d-wave symmetry of the order 
parameter. As far as modulus fluctuations in the particle-hole channel are concerned, in the d-wave case this would 
reflect in a Hartree-Fock-like correction to the band dispersion. Since this makes the analytical treatment not viable, 
both for the coherence length and for SD itself we do not address this issue in the present work. However, at generic 
flUing the physics is not expected to be different from the s-wave case, and the strong-coupling mapping of the effective 
action in the XY model should be preserved. 



VII. CONCLUSIONS 



In this paper we analyzed in detail the issue of the collective-mode description of a superconductor in the quantum 
regime. In particular, we addressed the issue from the point of view of the phase-fluctuation correction to the supcrfluid 
density, and we extended the previously known results for continuum microscopic models to the more realistic lattice 
case. 

Usually, the evaluation of the phase-fluctuation contribution to the superfluid-density depletion is addressed within 
the quantum generalizations of the classical XY model, which however fails in providing a complete description of 
dynamical effects. Specifically, in the continuum system the finite correction to D down to T = derived within the 
XY model implies the explicit violation of Galilean invariance. As we showed, the inadequacy of the quantum XY 
model can be understood as one derives the phase-only action starting from the microscopic model. This approach 
allows us to treat on the same footing both static and dynamic interaction terms in the phase field: besides the 
classical (static) anharmonic terms (V^)"*, one has third- and fourth-order quantum (dynamic) interaction terms 
which contain the time derivative of 6. These quantum terms are absent in the quantum XY model, where the 
dynamics only appears at the Gaussian level, and induce a correction to D which cancels exactly, in the continuum 
case, the contribution due to the classical interaction, restoring the equality ps — p at T — 0. We point out that this 
evident failure of the quantum XY in the continuum case naturally poses the question of the quantitative validity of 
the quantum XY model estimates even in the case of a lattice microscopic model. 

The detailed description of the main steps required to derive the phase-only action was motivated by the fact 
that this formalism also provides a more general framework for the description of different microscopic systems. As 
a consequence, we were able to investigate the phase mode both in the presence and in the absence of long-range 
Coulomb interactions, and for both s-wave and d-wave symmetry of the order parameter. As far as the lattice case 
is concerned, we also showed that the phase-only action acquires a more complex structure due to the fact that the 
minimal substitution couples the fermions to the electromagnetic field at all the orders. Moreover, since the stiffness 
at T = is no longer related to the particle density but rather proportional to the average kinetic energy (for nearest- 
neighbor hopping), no conservation law is violated by finding that classical and quantum interaction terms lead now 
to a finite one-loop correction 6D. When compared with the result of the quantum XY model, one also sees that 
6D is of order l/(fci?^o)^ with respect to 6Dxy- A similar result also holds in the classical regime, where we found 
that the phase-fluctuation correction to D is smaller (in the weak-to-intermediate-coupling regime) than within the 
classical XY model by the same factor ^ l/(fcFCo)^: in both the continuum and the lattice model. The reduction of 
the phase-fluctuation effects for kp^Q 3> 1 is made even more pronounced in the continuum case by the inclusion of 
long-range Coulomb forces. 

We also devoted particular emphasis to the discussion of the strong-coupling regime. As the interaction strength 
increases the fermionic model evolves toward a bosonic one, where the Cooper pairs act as individual bosons. Since in 
the superfluid the order parameter is connected to density, one expects that the description of the strongly-interacting 
electronic system requires the inclusion of both density and modulus fluctuations. We explicitly considered the 
negative-C/ Hubbard (s-wave) model. This model and its extensions are widely considered in the literature, as candi- 
dates to capture the main features of the intermediate-to-strong-coupling superconductors. This issue has attracted 
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renewed interest due to its possible relevance to the description of the underdoped regime of the superconducting 
cuprates. 

The inclusion of density fluctuations within the negative-C/ Hubbard model is performed by RPA resummation in 
the particle-hole channel of the fermionic bubbles which appear as coefhcients of the phase-only action. As far as the 
modulus fluctuations are concerned, their inclusion in the fermionic bubbles is directly derived by the definition of the 
susceptibilities in terms of the physical quantities, although the formal treatment of modulus fluctuation at Gaussian 
level is possible. The discussion of the result strongly involves the knowledge of the behavior of the coherence length 
^0 at strong coupling as a function of the fermion density, since it sets the spatial cut-off for phase fluctuations. This 
issue was addressed elsewhereji^ and allows us to conclude that at strong coupling, both in the presence and in the 
absence of long-range Coulomb forces, the fermionic system maps onto the quantum XY model, except in the very 
low-density quantum regime. 

Even though our attention was devoted to the correction to the superfluid density, our results for the structure of the 
phase-only action also suggest the possible outcomes of the analysis of phase-fluctuation effects on other quantities. 
The connection between the coefhcients of the phase-only action and the various physical response functions (in 
specific regimes) allows us to relate the correction to the Gaussian phase-only action to the phase-fluctuation induced 
correction to the corresponding physical quantity. This has been done, e.g., for the optical conductivity, as discussed in 
Ref. 30] in connection to the physics of high- Tc superconductors. A second possibility concerns the phase-fluctuation 
effect on the thermal conductivity, which is presently under investigation^. 
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APPENDIX A: DEDUCTION OF THE PHASE-ONLY ACTION BY MEANS OF FERMIONIC 

VARIABLE INTEGRATIONS 

In this appendix we report some details on the derivation of the phase-only action. As a first point we explicitly 
perform here the Hubbard- Stratonovich transformation on the microscopic action of Eq. |(SJ, introducing the auxiliary 
field A(q, r). In this way we get the decoupling of the interacting part Si of the action © in the particle-particle 
channel, which now reads 



S'/(co-,c+, A, A*) = E ~ V ^ A(q,rJm)w(k)c+^a^(w„)c+^^a^(u;„ - O™) 

^ \ffi J2 [A(q,r2„)]*u;(k)c_k+a|(u;„)ck+a|(w„ - (Al) 

Here A(q, fim) is the Fourier transform of the auxiliary field, q is the momentum and Qm — inm/ f3 is the Matsubara 
frequency. 

A second issue to be discussed is the effect of the gauge transformation 

c.(r,T)^4(r,r) = c.(r, T)e^'^(-^'-)/2^ 
c+(r,r)^c+:(r,r) = c+(r, T)e-^''(■^^-)/^ (A2) 

on the fermionic field appearing 5*0 and Sj. We start by rewriting Eq. (jAl|) in real space 

S'/(c^,c+,A,A*) === J y" l^i5llll!dRdT - J y"A(R,T)u;(r)c+(R+ ^,t)c+(R- |,r)dRdrdT 

A*(R,T)u;(r)q(R- |,T)cT(R+|,T)dRdrdT, (A3) 
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where w{r) = it;(k)e''''" is the Fourier transform of z«(k). After the transformation (jA2|) . we obtain, e.g., for the 
first term linear in A in Eq. (|A3() 



/ r|A(R,r)|e*^(^'")ii;(r)c+( 
Jn Jq 



R + - , r)c+(R - -,T)d'RdrdT - 

2 ^ 2 



I A(R, r)|e^[''(^^-) ^ ^1 u;(r)c+T (R + ^,t)c+^{R - r) 



dRrfrdr. 



(A4) 



ff the interaction is local in real space, i.e., w(r) = S{r), as in the (isotropic) s-wave case, the phase 9 disappears 

from the exponential in Eq. ljA4|l . In the c?-wave case instead, one can expand e'L''^"-^-! t' t^I ~ 1 - 

ir^r^Va V/3 6'(R) + • • •■ The residual dependence on the phase affects the form of the phase-fluctuation propagator by 
introducing terms which are irrelevant in the hydrodynamical limit. As a consequence, in the following we completely 
discard the dependence on the phase field in Si, assuming that after the gauge transformation (|A2|) Sj has the form 
(|X3|) with both A(x) and A* (a;) substituted by |A(x)|. 
As far as Sq in Eq. jSJ is concerned, the effect of the gauge transformation can be better expressed by first rewriting 



So = dx (f>'^{x) 



dr 







a. + 1- + M 



(f){x), 



where we introduced the standard Nambu spinor (/)+(r, r) — |^c^(r, r) ci{r,T)j. After the transformation (|A2|) . Sq 
gets modified into 



•^o = 5*0 + dx (f)'^ (x) 



idr9{x) {Ve{x)f 



8m 



T3 



Am 



Ve{x)- V 



To \ (kix), 



(A5) 



where the operator V= (V — V) acts on the fermionic variables. 

Introducing the Nambu notation also for the interacting part of the action, after the gauge transformation we find 



k'k ij q 



u 



(A6) 



where we put k = (k, cj„) and q — (q, ilm), and the matrix A^^,/. can be deduced from Eqs. (|Aip and HA5p . As we 
observed in Sec. Ill Al the fermionic variables appear in the action HA6|I only in the Gaussian term, whose coefficient 
^fc'fc contains the collective variables |A| and 9. By integrating out the fcrmionS(2S 



V^Vcp+e-^^'^^'^^'^' ^^^^^ 



DetA!^,^, 

k k 



and using the well-known identity: InDetA^"',^ = Trln A^^,^, where the trace Tr acts in the Nambu space and on the 



four-momenta fc, A:', we finally get the effective action for the bosonic degrees of freedom 

|A(g)P 



5e//(|A|,0)=5] 



u 



Trln^y,,. 

k k 



(A7) 



The second term of Eq. (|A7p acquires a more readable form by separating the part of the matrix Ak^k' with the 
explicit structure of 5kk' from the rest. This leads to 



TrlnA'^,, =Trln 

k k 



= Trln^o"^ + Trln 



= Trlngo-'-E^i//' 



N 



where 



TV times 



^eff - ^ 



^oS X • • • X GoJ: 
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and 




1 T' 

- ] e(k -k) V(k' - k - s) • se{k -k- s)e{s) 



The BCS matrix Qq is 



where 



|k|)(|k'| 



^0 



S 

|k|) 6{k - k) 



T3 



To 



k' 



|A|(A: -k) 



C/o(k,t^„) = 



•?^o(k,t^„) 



(k)VS|A|(0) 



El 



(A8) 



(A9) 



(AlO) 



and i?k = VCk + ■i«2(k)AQ. The value of the gap ampUtude A p = ^/T/fl\A\{0) in Eqs. (|A10|I is determined by the 
solution of the saddle-point equation for the effective action HA7|I . {dS/d\A\{q ~ 0)) = 0, and coincides with the 
mean-field value of the superconducting order parameter. The saddle point equation imposes a constraint only on 
the modulus of the order parameter, since a constant value of the Goldstone field 9 does not affect the action. As a 
consequence, Go and are exactly the BCS normal and anomalous Green functions for the Hamiltonian — 
By separating the field |A|(g) = |A|(0) -|-(5|A|(g) in its constant and fluctuating parts, the effective action HA7|) can be 
expanded in terms of phase and modulus fluctuations with respect to the BCS state. Observe that since one separates 
modulus and phase |A|(g) is the Fourier transform of the real field |A|(a;), and consequently [jAKq)]* = |A|(— g) and 
|A((;)p = J2q l|A|(9)P — 1^1(9)1^1(^9)1 where A(g) is the Fourier transform of the complex field A(a;). 
When, at the end, the gaussian expansion in the phase fluctuations around the saddle-point solution is performed, 
the Gaussian effective action reported in Eq. Ijf 0|l is obtained. 



neutral 



where the explicit expressions of coefficient A are the following 



with 



and 



App(9) 



T 



-^tr[go(fc)f3^o(fc + '7)T3] 



n 

k 

- > — tr 

k 



t?o(k+ ^,a;„)fo^o(k- ^,c^„ - ri„j)f3 



-^<5a. + A^^(,), 



A??(9) 



Pmf 



T k"^ k*" 
— > tr 

k 



^o(k+ - ,w„)fo^o(k - 



^m)To 



T 



51 tr [^o(fc)f3 



1 ^lltanh 



El, 

k 



(All) 
(A12) 
(A13) 

(A14) 
(A15) 



is the saddle-point fermion density. 

In the presence of long-range Coulomb interaction, the Hubbard- Stratonovich transformation of the microscopic 
action 10 reads 



Si (c,,c+,A,A*,ph5) 
- A*(g)y|^u;(kK 



~ V ^ ^ ^(k)<+.T('^»)^-k+ai(^" 



u(^n)ck+aT(^« 



■iPHs{q)]J cj_q,^(t^n - f^m)ck,<T(w„) 
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where the first, the third and the fourth terms were already present in Eq. ljAl|l . The last term describes the 
interaction of the fermions with the Phs- This term is not affected by the gauge transformation l|A2|l . since it is of 
the form c+c. Integrating out the fermions the effective action for the collective variables I) A7(l is now 

5.„(|A|,«,™,,) = E + ^} - T'>"4v 

where the matrix elements J^^,^. include now the density field phs, ^k'k — * ^k'k + iyT/VLpHs{k — k )t^. As a 
consequence, the density contributes to both the diagonal part of J^^^., whence 

^o"' ^ ^o"' - *y|pH5(0)f3, (A17) 

and to the self-energy matrix S — > E + ipHs{x)T^- Again, we separate phsIq) — PHsiO) + SpHsiq)- The value of 
Phs{0) is determined by the solution of the saddle-point equation for the action (|A16|I . i.e. 

''"'^'^ -^ j:iSo{k) + go{-k)]^PMF 



y(o) n 



k 



where pmf is the particle density evaluated at mean-field level according to Eq. ljA15(l . Therefore, the correction 
(|A17(I to the Green function is the standard Hartree shift of the chemical potential /j, /i + V{0)pmf -^^ ■ 

In the charged case, as anticipated in Sec. II B, in order to get the phase-only effective action at the Gaussian level, 
it is necessary at first to evaluate the Gaussian expansion of 5*6// both in the phase and the density fluctuations 



Oi~q) 

^PHsiq) 



where 



B{q) = 



nl^Appiq) - qaClbAfjiq) + 2inmClaApj ~2ciaKj{q) + 2in^App{q) 
2q,A;j(q) - 2in„,App{q) 4 (pJ^ + App(q)') 



with the same definitions (|All|) - (|A13p of the coefficients A, and, in a second time, to perform the Gaussian integral 
over the field Spns- 

As a final result of this procedure, the following expression Hl()|l for the effective action is obtained 

q 

where the coefficients C are given by 

i^pp{q) = 



Appjq) 

^'^'"^ l + Vi<l)Appiqy 



KM 



^^^^'^^ l + V{ci)App{q) 

ab / \ A ab 



a;,,(?)aj,,(?) 



APPENDIX B: EFFECTIVE ACTION AND GAUGE INVARIANCE 

In Sec. Ill Al we observed that the BCS bubbles A{q) which appear as coefficients of the Gaussian effective action 110(1 
break, as it is well knowni^, the gauge invariance of the theory. In this appendix we derive the general relationship 
between the coefficents A{q) and the corresponding gauge-invariant electromagnetic response functions K{q), which 
control the physically accessible quantities (see also Ref. (511131]'). 
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We first introduce the electromagnetic potential = (0, A) into the Gaussian effective action (|10|l via the minimal 
substitution {86 /dt) {89 /dt) + 2etp, Vfl -> V6I - (2e/c)A, which gives 



e2 



Saie) - Saie^A) = Scie) + - Y,A^{q)K^^{q)A,{-q) + Y,^e'[q^t^^^{q)AMO{-<l) - q^.K>^''{-q)A,{-q)e{q% 



g 



where the four-dimensional notation for A^'^ was introduced, so that App — A°", A'J^j — and A'fj — A°^. We can 
now integrate out the 9 field to obtain the partition function Z[A^] which allows one to define the electromagnetic 
response functions of the system as 



K^^{q)- 



5A^{q)5A^{-q) 



= A-(,) + ^^^:^(%^^. (Bl) 



This equation provides the relation between the coefficients A^'' of Eq. 11U|) and the physical response functions K^^'^ . 
Now it can be easily seen that the difference between the A and the K functions (i.e., the second term in the right- 
hand side of Eq. (|B1|) ) has a purely longitudinal structure. In other words phase fluctuations affect only longitudinal 
correlation functions^^. As a consequence, we can safely identify the transverse part of the Gaussian bubbles (A"'") 
with the transverse part of the physical response functions [K'^)- 

This fact turns out to be crucial in interpreting the Gaussian action (|10l) in the hydrodynamic limit, which is 
the physically interesting case while evaluating the low-temperature behavior of the collective mode. Indeed, by 
calculating the static limit = 0,q — > of the coefficients (|A11|) - (|A13|I we observe that: 

• In the Gaussian action (|10|) the longitudinal part of Ajj appears, which, as we just showed, is not the physical 
longitudinal correlation function. However, in the q ^ limit, the transverse Ajj and longitudinal Ajj part 
of Ajj do coincide, leading to the identification of the coefficient of with the stiffness, defined through the 
transverse physical function (see, e.g., Ref. 

lim A^j(q, a„ = 0) = lim Ajj(q, a„ = 0) = lim Kjj{q, n„r = 0) = -D{T). 

q— >0 q— »0 q-^0 

• The coefficient A°j(q) vanishes in the static limit, as it can be seen by considering its definition lA12p . 

• From Eq. (jBip it turns out that the coefficient App coincides in the static limit with the physical density-density 
correlation function Kpp, whose static limit is by definition the (bare) compressibility Ka{T) of the system, 
limq^o App(q, fi,„ = 0) = limq^o Kpp{q, f2„ = 0) = ko{T). 

As a consequence of the above results, the Gaussian phase-only action in the hydrodynamic regime as exactly the 
expression presented in Eq. (|10|l . 
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